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Abstract 

Long, elastic filaments cross-linked and deformed by active molecular motors occur in various natural settings. The 
overall macroscopic mechanical response of such a composite network depends on the coupling between the active 
and the passive properties of the underlying constituents and nonlocal interactions between different parts of the 
composite. In a simple one dimensional system, using a mean field model, it has been shown that the combination 
of motor activity and finite filament extensibility yields a persistence length scale over which strain decays. Here 
we study a similar system, in the complementary limit of strong noise and moderate extensibility, using Brownian 
multi-particle collision dynamics-based numerical simulations that includes the coupling between motor kinetics 
and local filament extensibility. While the numerical model shows deviations from the mean field predictions due to 
the presence of strong active noise caused by the variations in individual motor activity, several qualitative features 
are still retained. Specifically, for fixed motor attachment and detachment rates, the decay is length is set by the 
ratio of the passive elasticity to the active shear resistance generated by attached motors. Our study generalizes 
the notion of persistence in passive thermal systems to actively driven systems with testable predictions. 
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Introduction 

The active strain generated by molecular motors moving 
on elastic filaments is the principal mechanism of force 
generation and motion in cell biology [BlllllSlillZlIl 
[5]. In natural and reconstituted versions of such active 
composites, the localized strains due to motors deform 
the filaments dynamically, due to the interplay between 
elasticity, geometric constraints, active motor forces and 
noise due to fluctuations in the ambient medium and 
in motor activity. These deformations in turn regu¬ 
late the activity of the motors themselves mils]. A 
structurally ordered example of such a system is the 
eukaryotic flagellum that is made of relatively stiff fila¬ 
ments (microtubules), motors (dynein) and passive elas¬ 
tic elements (nexins) which together oscillate with well 
defined wavelengths and frequencies [HIIIIIIIIIIIIIII. 
While a variety of models of varying degree of com¬ 
plexity [HElEllzlElIH] are consistent with observations 
of increasing wavelength with increasing flagella length, 
experimental evidences indicate that this relationship 
is altered for very long lengths [13 Hi [m ng and 


that the wavelengths attained in nature are self-limiting 
Oil] even as the flagella themselves range from tens 
of microns to nearly a centimeter. This strongly sug¬ 
gests that over large lengths, mechanical information 
transmission degrades substantially. 

In a passive context (3 , it has been recently demon¬ 
strated that when bundles of filaments are forced to 
bend, the shearing forces between them are mediated by 
extensibility, leading to a characteristic scale over which 
mechanical information is transmitted. In an active sys¬ 
tem such as an array of soft filaments driven by molec¬ 
ular motors, a possible consequence of this finite decay 
length of strain is that the weak extensibility of the fila¬ 
ments can limit the range of mechanical signal transmis¬ 
sion between molecular motors and restrict the length 
scale over which motor coordination can occur. This 
will significantly affect the motor activity in filaments 
of lengths much larger than the strain decay length and 
could bring in spatial inhomogeneities in both strain 
and activity. As an example, in Fig. 1(a), we show two 
distinct motor patches labeled I and II cross-linking a 
pair of thin filaments that may bend and shear. For 
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Figure 1: (a) An active, motor-filament composite illus¬ 
trating two distinct motor aggregates, groups of motors 
(I) and (II). Of interest is the distance between two 
motor groups beyond which they act as if isolated, (b) 
Schematic and geometry of the passive filament com¬ 
posite with length effective thickness D + 2w. The 
arc-length is measured along the tangent t. (c) A sketch 
of the coarse-grained active gel (top) corresponding to 
the mean field approximation of the microscopic fila¬ 
ment / motor system (bottom) illustrating the motors 
(blue) and cross-linking passive linkers (red). 


inextensible filaments, mechanical activity by group (I) 
and shear (sliding) induced by this patch is transmitted 
by attached motors in (II) over arbitrarily large inter¬ 
aggregate distances. For extensible filaments, there is an 
interplay between shear / slide and filament elongation 
along its contour length and leads to a decrease in the 
ability of group (I) to mechanically link with group (II). 

Using a mean field approach valid in the noise-less 
limit [To], we recently showed that the decay length of 
localized steady or oscillatory strains is determined by 
a combination of filament elasticity, passive shear resis¬ 
tance imposed by cross-linkers and the active viscoelas¬ 
tic properties of the motor aggregates. This analysis 
was performed by coarse-graining the stochastic nature 
of motor-filament interactions, and ignoring variations 
in motors kinetics due to the finite number of motors. 
However, even if thermal noise can be ignored (such as 
when motors are rigidly fixed at the base and there¬ 
fore do not diffuse) there is still an effective motor noise 
that arises due to discreteness of binding and unbind¬ 
ing events and fluctuations in the attachment times due 
to the finite number of motors. Moreover, there is a 


complete two way coupling between the motor activity 
and filament dynamics that arises as the motor activ¬ 
ity causes filament strain, influences the collective motor 
dynamics. All these effects are not captured in the mean 
field model. 

Here, we present a microscopic numerical model, in 
which we incorporate the effect noise arising from vari¬ 
ations in motor kinetics, and implement the two-way 
coupling between motors and filament to study its influ¬ 
ence on the persistence of active strain by comparing 
the results with the mean field predictions. Predictions 
of the microscopic model agree qualitatively with the 
mean-field results, despite the effect of strong active 
noise. However, motor noise is shown to result in richer 
dynamical features including modification of the motor 
duty ratio and localized regions of coherent oscillations 
due to the interplay between the local active shear 
resistance and filament elasticity. 


Persistence of strain in a passive composite 

To understand the coupling between shear, bending and 
extension in a passive setting, we start with a mini¬ 
mal model of two elastic filaments of unstrained length 
I, lateral width b and thickness w i held at a dis¬ 
tance D <^ £ apart by a series of passive, linearly 
elastic springs with areal density pN and stiffness An 
(F ig. Sb) -(c)). The passive springs are compliant in 
shear along the axial direction but prevent any change 
in the distance between the filaments. 

At length scales large compared to the spacing 
between these springs, the composite acts as a filament 
of thickness D, with an effective shear modulus that 
depends on the spring stiffness, the density of springs, 
and the inter-spring spacing when b ^ max(iy,ZI). By 
considering a sliding displacement of the top filament 
relative to the bottom due to a shear acting on its sur¬ 
face of area b£, with a stress as, the net force is agbl. 
Further, if the areal density of springs is pni the shear is 
balanced by p^bi springs, all acting at the same time, 
and each contributing to shear resistance with spring 
constant An- Assuming a shear (sliding) displacement 
of magnitude A, force balance 

G*A 

A kyi b i as b£ ^ ^ b £ 

yields an effective shear modulus from these passive 
springs given by G* ~ HAnPn- For a strip of lateral 
width w, we have G*w/D ~ wAnPn - consistent with 
previous analyses of a passive cross-linked railway track 
model [ZIIS]. 

With this estimate in mind, we choose a local coor¬ 
dinate system characterized by an arc-length variable s 
so that the filament is in the range s G (0,£). For two 
points ri (on filament 1) and r 2 (on filament 2) (Fig. 
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Bb)) which face each other when the filaments are in 
the undeformed state, we define the sliding deformation 
A(s) as difference between ri and r 2 , relative to its ini¬ 
tial difference before deformation. The angle made by 
the tangent to the centreline, 6{s) is then related to the 
slide A(s) and the displacement u{s) by geometry via 
the relations dri/ds = t, dt/ds = Cn and dn/ds = —Ct, 
where n and t are the normal and tangent vectors, and 
C = dt/ds-n « 9s is the local curvature. In the deformed 
state, the position r 2 relative to ri can be written as, 

r 2 « ri -f Dn-fut. (1) 

Here we have assumed {D + w) £. The axial strain 
along the filament is Us where d{.)/ds = (.)s. Using the 
relations A^ = i:)(|(ri)s| - |(r 2 )s|) and (ri)^ = t, we 
obtain, 

(r 2 )s ^t — DCt +Us t + uCn. (2) 

For small deformations we can ignore quadratic/cross 
terms in the deformation to get 

|(ri)s| = 1, |(r2)s| « 1 - + Us (3) 

Substituting ^ into ([^ we get 

As = Dds-Us. ( 4 ) 

Integrating Q from s = 0 to s using the boundary 
condition 0(0) = 0 yields A(s) = D9{s) — u{s). 

The shape of the filament is then obtained by mini¬ 
mizing the energy Et of the composite due to bending, 
stretching and shear, with 

Bpasd^ + KpasUg + GpasA^ 

2 

Here Upas Ew^, fsTpas ~ Ew and Gpas ~ Pn^n are 
the passive bending, stretching and shear moduli, with 
E is the Young’s modulus of the material. The Euler- 
Lagrange equations obtained by minimizing the func¬ 
tional are given by 6 Ex/S A = 0 and 5Ex/59 = 0 
and lead to equilibrium static solutions. 

We consider solutions to Q in two complementary 
limits. When the bending stiffness Upas —>■ co, the com¬ 
posite deforms due to shear and extension alone and 
A = —u, so that the minimization of provides 

A,, - (Gpas/i^pas) A = 0, (6) 

which yields the relaxation length scale for pure exten¬ 
sion £\ = -f^pas/Gpas Ew/pjqkn- For inextensible £1- 
aments, i^pas oo, and A = DO and the minimization 
procedure provides the relation 0ss— (^^Gpas/-Bpas) 0 = 
0, thus yielding the persistence scale for pure bending 
£% = Hpas/T>^Gpas ~ Uw^/H^/cnPn- 


ds. (5) 


Ex — b 


Allowing for sliding deformations at s = 0 and 
keeping bending, extension as well as sliding terms, we 
repeat the minimization procedure to obtain the two 
coupled equations, 


Bpas^ss + DKpas{Ddss — Agg) — 0 (7) 

Kpas{Ass — DOss) — GpasA = 0 (8) 


In this general case, the second equation can rewritten 
as Ass — £ 7 '^ A = 0 and indicates that variations in the 
sliding displacement A are associated with a persistence 
(decay) length 




£b£e 

\J£\e £\ 


w/D 


{KpaslGpa.Y . (9) 


Similar length scales appear in a variety of soft systems 
where bending, shear and extension are coupled, both 
in microscopic and macroscopic contexts mm- 

To put these length scales in perspective, we look at 
typical parameters for a flagellar axoneme. We take w 
and D as the radius and the spacing of microtubules, 
E as the Young’s modulus of microtubules, as the 
stiffness of passive (nexin) links and pN as the density 
of these links. Using w ^ 20nm, D ^ 40nm [3], E ^ 1.2 
GPa [21 Ej, 16 — 100 pN pm“^ [23] and wp^ ~ 
10® — 10^ m“^ [T71122] I we estimate £e ~ 200 — 500 
/im while ^ 80 — 200 pm. This estimate is however 
with just the passive cross linkers contributing to the 
shear stiffness and therefore values should be treated as 
an upper limit. Motor activity will naturally influence 
this scale as attached motors contribute to the instan¬ 
taneous shear resistance. Assuming that all motors are 
bound, with linear density 0(10®) m“^ and effective 
spring stiffness of 10“® N/m, we obtain £e ~ 5—10 
pm. Since sperm flagella are often much longer than 
this scale [4], the role played by activity in modulating 
persistence lengths can be important. 


Mean-field model for an active composite 

For the case we treat here, corresponding to the case 
without bending wherein the slide A(s) = —u(s) 
according to Q, motor activity modifies the decay 
length £e in the noiseless limit by changing the effec¬ 
tive shear resistance in a manner controlled by motor 
kinetics. We being by summarizing the results of the 
mean held analysis |10j , valid in the noise-less limit with 
the focus on systems with purely extensional modes of 
deformation. While the original analysis m allows for 
a more general form for the motor-filament interactions, 
here we consider a special case of the continuum theory 
consistent with our simulations as shown in FigQb). 

Here, two cross-linked filaments are held together 
by a combination of active elastic motors and by pas¬ 
sive elastic spring-like linkers. The active motors act as 
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linear elastic springs of stiffness k^, and can either be 
attached to or detached from the filament. When the 
inter-link spacing and the inter-motor spacing is much 
smaller than the filament length, the motor-filament 
assembly can be represented by a continuum descrip¬ 
tion. In this case, for an inertialess filament, the sum of 
forces on the filament arising from filament extensibil¬ 
ity, springs and motors must vanish. If F^n is the active 
force due to a single attached motor and pm the motor 
density, this yields 

i^pas'^s)s Gpasri 4“ Pm^m — 0; (^0) 

which differs from previous system described by 
because of the active force density term PmFm- 

To estimate the active force density Fm that con¬ 
tributes to an active shear resistance, we consider a 
minimal model for the motor-filament interaction. In 
this, active motors attach to the filament with zero 
extension (no pre-strain) and once attached walk rel¬ 
ative to the filament with a velocity Vm- The motor 
velocity Vm is linearly related to strain through a linear 
relation, Vm = yt = Ut + ?;o(l - k^y/Fs), where Ut is 
the velocity of the filament relative to the motor, Fg is 
the motor stall force and Vg is the zero load velocity of 
the motor. Since one end of the motor is permanently 
attached to a rigid substrate, motor movement results 
in a motor extension y and thereby causes a resisting 
force kraV- 

Motor kinetics may be described by a set of popula¬ 
tion balances relating the attached and detached prob¬ 
ability densities to the corresponding fluxes. Ensemble 
averaging the terms in these relations provide dynami¬ 
cal equations for the fraction of attached motors, N{s,t) 
and the mean motor extension, Y{s,t) (see [TU] for 
detailed derivation). The evolution equation for N is 
then given by 


Nt = -N)- ccoffiv. 


( 11 ) 


where a;°„ is the mean attachment rate of motors, and 
0 is the load dependent detachment rate, defined as 
0 = 0°F{£,Y), where £ = k^S'^/ksT with being 
the extension corresponding to motor detachment, and 
the scaled mean extension of the attached motors, Y = 
{y)/^m- Since the motor detachment is strain depen¬ 
dent, the evolution of N is non-linearly coupled to the 


mean motor stretch Y, as implied in (11). Population 


balance relations then lead to an analogous evolution 
equation for Y that is given by 


kt — ^ -I- AiIjjIs {A 2 - 1^) + ^ ). (12) 


In (12), the first term on the right hand side is the 


stretch due to the passive motion of the attached motor. 


the second term gives the motor velocity relative to the 
filament, and the third term corresponds to the rate 
at which the mean strain changes due to the kinetics 
of motor attachment and arises from the difference in 
extension of attaching and detaching motors. The two 


dimensionless parameters in (12) Ai = Vgkm/{0°Fs) 


and A 2 = Fs/{km5m) relate the microscopic motor 
kinetics to the population averaged dynamics. Once we 


have N and F, the motor force Fm in (10) is written as 


p — i- A AfY 

± m. — '^mOrri 1» J • 


(13) 


Equations (10)-(13) provides a closed system for the the 


dynamics of filament displacement u, attached motor 
density N and motor extension Y in terms of dimen¬ 
sionless parameters Ax, A 2 and 4' = 0/won- These 
equations have a homogeneous stationary solution given 
by A^o = (1 + 'I'-^o)”^ and Yq = AiA 2 iAi + Fo)~^, and 
a natural question is that of their stability to variations 
in the parameters. 

Before we describe the global distribution of the 
strain field and the persistence of activity, we first exam¬ 
ine the dynamics of a small fragment of the filament- 
motor composite of length Ig ^ (ATpas/Gpas)®- In 
this limit, we can ignore filament extensibility and the 
motion of this segment relative to its neighbors can 
be mapped to that of a homogeneous population of 
motors acting on a rigid segment and working against 
an external spring with effective stiffness Kg oc GpasQ 
Eor such a fragment, when Kg > 0, a linear stability 
analysis of equations ([n])-(i^ reveals that the station¬ 
ary states JVq and Yq become linearly unstable and a 
stable oscillatory state emanate via supercritical Hopf- 
Poincare bifurcations [TU]. In the stable oscillatory state, 
the power input into the system due to activity bal¬ 
ances the energy dissipated by motor viscosity. The 
effective motor friction characterizing this process is 
Vact ~ PmNgkm{d{y)/dut)ut=o where PmNg is the total 
number of attached motors. As 4/ increases, the term 
{d{y)/dut)ut=o typically decreases, as does Nq thus pre¬ 
dicting an increase in the frequency of the emergent 
oscillations uic- 

To connect these localized solutions to the global 
extensional field in order to determine decay lengths of 
imposed extensions, we consider an active composite fil¬ 
ament held fixed at one end (s = 0) while free at the 
other (s = t) which is forced harmonically. We assume 
that u{0,t) = 0 and u{£,t) = ug{tj + euReal[e“*], 
where eu/uo{£) 1 characterizes the deviation from 
the stationary extensional state Uo{s) which satisfies the 
equation 


Aipas('ao)ss GpasriO 4” GactG — 0, 


(14) 


^Strictly speaking, there is also an active contribution to this 
effective spring constant that comes from attached motors in 
neighboring filaments - but this just rescales Ks and the mapping 
still holds. If pN = 0 then the passive part of Ks oc i^pas- 


Biophysical Journal 00(00) 1—10 











Persistence BJ Submit — 2015/1/13 — 1:36 — page 5 — #5 


5 


where G^ct = Pmkmi the active analogue of passive 

-4^- With 


shear modulus and C = Sm^oYo = ^ ■ 

boundary conditions uo(0) = 0 and (ito)s(V) = 0, at 


leading order, i.e. 0(e°) (141 provides, 


uo(s) = /3 ^1 — a e^E — (i — q-) e ^ 


(15) 


Here, we note that the decay length Ie = 'v/'^pas/Gpas, 
a = e-"^(l + and /? = 

does not depend on the active shear resistance because 
motors sense only strain rates and not the actual strain. 

At the next order, i.e. 0(e) we find the relations 
that determine the linearized response of the system to 
externally imposed small amplitude perturbations. The 
magnitude of the time dependent part of the extension, 
u{s) due to the active force from the motors PmFm can 


be found by substituting PmFm = —GactX^ in ( |T0| and 
yields 


AT, 


pas ^ss 


(G^pas “1“ GactX) ^ 


(16) 


where x(w) = —SmiYoT 2 + NqTi) characterizes 
the complex linearized response of the system 
with Ti = - {iu+ u}°^'liFo+u}°gAi) - 

' iui+ui °^+ui°lfT o ' 
iiT'No 


(e) (<.. (H) (= 


0 ) 


-1 


and To = 


off^O^O 

equation of the form u 
expression for the effective persistence length. 


Seeking solutions to this 
.j gs/A yieifjs following 


Xe = Real[A] = V2 ( \Jp\ + P 2 + Vi 


(17) 


s/Aipas) “t“ (G^act/ATpas) 7^ and <^2 — 


(Gact/ATpas) T with TZ and I being the real and imagi- 


where Lpi = (Gpas/ 

)I 

nary part of x respectively. We note that in the absence 
of activity, i.e. Gact = 0, we recover the expression 
X = £e- 

In summary, the mean-field analysis predicts that 
the relative importance of the active (motor) to passive 
(linker) elasticity is controlled by Gact/Gpas, whereas 
in the absence of passive linkers the persistence length 
is controlled solely by Gact/ATpas- Further, (171 implies 


that the persistence length is influenced by the imposed 
frequency w, though the frequency dependent complex 
response term, x that controls the competition between 
intrinsic motor duty ratios and the extrinsic imposed 
frequency ui in determining the filament extension. 


Numerical model for the motor-filament composite 

However, the mean field theory does not take into 
account the effect of local motor noise due to dis¬ 
crete attachment and detachment events and the finite 
number of motors. Also, there is no intrinsic coupling 


between filament extensibility and the mean attachment 
time of motors - as evidenced by the fact that Nq and 
Yq are independent of uq. We now extend this minimal 
theory and investigate the role of increasing noise and 
fluctuations on the decay of both steady and oscillatory 
localized extensions. 

In our model, the continuous filament is made up of 
N discrete spherical monomers of diameter a located at 
Tj, (i = 1, ...jN), connected by an elastic potential 

Af-l 

Y H 

Bending stiffness is implemented via a three-body bend¬ 
ing potential 

N-l 

wb = ^ y] (ti+i - ti)^ (19) 

i=2 

where tj = (ri — ri_i)/|ri — ri_i|. Here the effective 
bending rigidity k penalizes angular changes from a 
local straight geometry. We fix the filament length to 
be a constant £ = 80 a, and also hx the spring constants 
Ki = K = 2 X 10'^ {UbT/u"^) in order to get a large fil¬ 
ament persistence length ip ~ 250 i. The bending and 
extensional stiffness are large enough so that the seg¬ 
ment may be treated as a rigid inextensible segment. 
For simplicity and to focus on the role of activity and 
noise, we did not incorporate permanent passive cross 
links in the model. With this simplification, shearing is 
resisted because of the temporary crosslinks (bridges) 
formed by attached motors. These bridges disappear as 
motors detach at sufficiently high loads. 

The active motors were modeled as linear springs 
with stiffness km and equilibrium length im- The motors 
attach with zero mean strain and walk along the fila¬ 
ment with discrete step size dm- As the motor heads 
move, the accompanying extension of the motor length 
leads to a force = —Hm{i — im) i that acts on 
the attached filament. Thus, motor activity results in 
localized active strain on the section of the filament 
to which they are attached. Consistent with previous 
experimental studies the load dependent stepping rate, 
i.e. the stepping velocity of the motor head - is chosen 
to have the simple form Vmi^ — |Fm|//max)- However, 
unlike in mean-field models where the attachment and 
detachment processes are coarse-grained into rates, in 
the simulations we prescribe microscopic probabilities: 
the load independent attachment probabilityponj is the 
probability of a detached motor to attach to the fila¬ 
ment, while the motor detachment probability Poff is 
defined by the piecewise function 

if |F,„| < /cut 
if |F,„| > /cut. 
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Figure 2: (a) A schematic diagram of the filament-motor system for a small fragment. The filament moves due to 
active forces imposed by attached motors - this motion is resisted by a anchoring spring of stiffness Kg. (b) The 
time dependent displacement, Aa{t)/dm for a fixed value for the ratio of connecting spring to motor stiffnesses, 
Ks/km = 0.2, and for zero load motor detachment/attachment probabilities, Pos/Pon = 1-0 (top), 2 (middle) 
and 10 (bottom), illustrating various dynamical regimes observed in simulations, (c) Phase diagram indicating the 
dynamical regimes as a function of the ratio Kg/km, and the ratio of zero load detachment-attachment probabili¬ 
ties (poff/Pon)-Tlae observed dynamical regimes are, steady extension (o), intermittent oscillations (♦), and steady 
oscillations (■). Dashed curves are a guide to the eye. (d) The oscillation frequency, Wc of a periodically oscillating 
filament with Kg/km = 0.5, as a function of Pos/Pon- 


This load dependence of the detachment probability of 
individual motors allows for the crucial two-way cou¬ 
pling between filament elasticity and the motor response 
and also sets the critical motor extension at which 
detachment occurs. We take /cut = /max in our sim¬ 
ulations, although in general they are not equal. Fur¬ 
thermore, the constant load independent detachment 
probability PoS in conjunction with pon yields a well 
defined equilibrium attachment of motors at zero load, 
such that, Pon(l — N) + Po^N = 1. However, we note 
that despite the added noise, our numerical model is 
an approximate one in that excluded volume interac¬ 
tions between the motors are not taken into account; 
these interactions will be significant only at large values 
of Pon when almost all motors are always attached, a 
situation that we do not consider here. 

The position and velocity of the monomer beads are 
updated using velocity-Verlet algorithm, in which the 
force acting on the monomer is calculated at every time 
step. The entire system is immersed in a viscous medium 
using the Brownian multiple particle collision dynamics 
(MPC) scheme which neglects the effects of inertia [Mj . 
According to this scheme, each monomer independently 
performs a stochastic collision with a phantom fluid 


particle with a momentum taken from the Maxwell- 
Boltzmann distribution with variance pkbT, where p is 
the fluid density. 

Using the numerical model, we first analyze the 
effect of noise in dynamical regimes of a short inexten- 
sible fragment - this corresponds to the dynamics of 
filaments with length ^s ^ (ATpas/Gpas)®- This allows 
us to identify the way noise modified the predictions of 
the mean field model, which predicts a sharp separa¬ 
tion between two states - a stable, stationary state and 
a stable, oscillatory state. We expect noise to make true 
boundary between these two regimes fuzzy. 

In Fig. [^a), we show a rigid filament anchored to 
a wall via a linear elastic spring of stiffness Kg that 
moves under the action of 800 motors. Monitoring the 
long time displacement of the free end as a function 
of the spring constant Kg and the probability ratio 
Pos/Pon yields results summarized in Fig. [^b)-(d). The 
interplay between the softness of the filament compos¬ 
ite Kg/km and the activity Pos/Pon yields three distinct 
dynamical regimes illustrated in Fig.j^b). For fixed Kg 
when Pos/pon 1; there is a large force on the fila¬ 
ment, and the free end of the filament attains a steady 
displacement with a well defined mean attached motor 
density, with fluctuations due to the noise imposed by 
discrete nature of motor binding and unbinding. For 
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Figure 3: (a) Schematic sketch of a composite filament made of 50 rigid segments, with s as the coarse-grained 
position measured along the fiament. (b)The mean steady state extension of material points as a function of their 
position from the clamped end (s = 0). The vertical bars correspond to the root mean square deviation and indicate 
the role of noise. The extension is less correlated for soft assemblies, when Kp^s/km = O.Ol(left), while it is highly 
correlated when Kp^s/km = 10. (c) Time averaged steady state extension of the segment at the free end, uo{i), 
as a function of zero load motor attachment/detachment probability, Pon/PoS, for various values of passive spring 
stiffness, Kp^^/km- The extension increases with Pon/Poff: and decreases with Kp^/km, consistent with analytical 
predictions. For large Pon/Poff the extension saturates at a maximum value, set by the finite length of the filament- 
motor composite, (d) The fraction of motors attached on a composite filament consists of 50 segments, as a function 
of their position from the clamped end for Kp^s/km = 10. The mean motor attachment is maximum close to the 
clamped end, where the local strain is relatively small. 


Poff/Pon ^ 1, the average number of motors attached to 
the filament is close to zero, and there is effectively no 
active force. However there is an intermediate range in 
Poff/Pon, where we observe first intermittent (Figj^b)- 
middle) and then regular (Figj^b)-bottom) hlament 


oscillations. While the regular oscillations are consis¬ 
tent with the mean held predictions, the intermittent 
oscillating state - where the hlament-motor composite 
randomly switches from a stationary to an oscillatory 
state and back, is a new dynamical state observed only 
in the presence of noise. A sweep of the parameter space 
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Figure 4: (a)-(c)The time averaged amplitude of imposed oscillations of the filaments as a function of scaled posi¬ 
tions, for various elasticity contrast (a) Jfpas/^m = 10“^, (b) 10“^ and (c) 10“^. An exponential decay length is 
indicated whose value depends on both motor kinetics as well as the elasticity contrast, (d) The decay lengths 
obtained by analyzing the amplitude of oscillations as a function of Kp^s/km for various Pon/PoS- 


yields the phase plot in Fig. [^c). The parameter range 
over which steady oscillations are seen is qualitatively 
consistent with theoretical predictions m- However, 
due to discrete motor noise, the boundary separating 
the steady and oscillatory states is no longer sharp - 
as expected, we observe a region of intermittent oscil¬ 
lations. Within the oscillatory regime, the frequency 
increases with increasing Pos/Pon (Fig. [^d)), which 
also agrees with the trend seen from the mean field 
predictions of the theoretical model. 

Next we study the role of motor noise on the exten¬ 
sibility of a composite filament that is made up of 50 
equally sized rigid segments, each linked to adjacent 
neighbors by linear springs of stiffness Ks (Fig. Sa)) 
that account for the weak extensibility of the compos¬ 
ite. Since no permanent cross links are present, pN = 0, 
and the effective passive stiffness Kg oc ATpas- We choose 
parameters for the filament elasticity and attachment 
and detachment probabilities such that the base state 
in the absence of imposed oscillations is steady (non- 
oscillatory) and stable and oscillate the free end at low 


frequency (w and small amplitude (Ai ^ 

Aa). For sufficiently low frequencies, the results are 
qualitatively independent of w and hence we will focus 
on results obtained for a frequency w = 10“^(drn,/'Cm)- 
Quantifying the mean extension of the filament as a 
function of arc length s of the extensible filament motor 
composite shown in Fig.j^b), we see that the maximum 
at the free end is influenced strongly by both Pos/Pon, 
and ATpas- However, the discrete nature of motor activity 
causes large fluctuations in mean extension, especially 
for very soft filaments with ATpas/^m ‘C 1 (Figj^ -left). 
Motor noise also results in uncorrelated spatial domains 
in extension for ATpas/^m = 0.01 while a greater degree 
of correlation is observed for ifpas/Am = 10 (Fig|^b)- 
right). Concomitantly, variations in the mean extension 
(shown as vertical bars) are also more correlated for 
stiffer filaments than for softer filaments. The noise 
decreases with increase in pon, as shown in Fig. |^b). 

Fig-i c) shows that the steady state extension of the 
free end uq {£) /a increases with Pon/PoS ■ To compare the 
numerical predictions for steady state extension mo(s) to 
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our mean field theory cni, we simplify by setting 
Pn = 0, yielding the solution 


uo(s) 


/ Gact / 

1 A1A2 , 

( 5m \ 

\KpaJ 

' Al + J-Q ' 

[l + ^ToJ 


s{2^ 


s) (20) 


from which one deduces that there is no decay length 
scale in the absence of permanent passive linkers. We 
see from the right panel in Fig. [^b) that this is needed 
the case. While the simple theory predicts a quadratic 
form for the extension consistent with Fig. I^b), we 
see that increasing the attachment probability yields a 
sharper gradient close to the clamped end and a flatter 
profile near the free end. Setting G^ct = Pmknn we find 
from ( |20[ ) that the extension is an increasing function 
of km/Kpas and is maximum at the free end, which are 
both consistent with our numerical results (Fig. ic)). 
The dependence on the activity 'k is also qualitatively 
captured as the slope increases with the ratio Pon/PoS- 
The extension at the free end UQ{i) obtained from 
(20) indicates that uo{£) decreases with Kpas/km while 
it increases with (pon/Poff) - both these trends are 
observed in the numerical simulations(Figj^d)). 

Since the motor detachment is strain dependent, the 
motor imposed extensions in turn regulate the motor 
kinetics - and specifically modify the mean attached 
time. In Fig. |^d) we plot the mean fraction of motors 
attached, (iVattached)/.^ 0 ) as a function of the distance 
from the anchored end, s for it'pas/^m = 10. When 
Pon/Poff ^ 0.1 the motor imposed extension of the fila¬ 
ment is weak due to less number of motors attached, and 
the distribution of attached motors along the compos¬ 
ite filament is approximately uniform. However for large 
values of Pon/Poff (Pon/PoS /t 1) the enhanced motor 
attachment imposes significant extension to the fila¬ 
ment, which is maximum at the free end. The nonuni¬ 
form extension of the composite leads to a nonuniform 
distribution of attached motors (Fig[^d)), where the 
mean motor attachment is maximum near the anchored 
end (s = 0) and decreases with the distance from that 
point, as the filament strain increases. This leads to a 
diminished motor activity near the free end. These fea¬ 
tures are not seen in the mean field predictions and are 
thus a direct consequence of the noise due to the discrete 
nature of motor binding and unbinding events. 

We finally calculate the length-scale corresponding 
to the decay of the amplitude of the imposed oscil¬ 
latory strain {u) along the composite, for a range of 
Kpa.s/km and Pon/Poff- For all values of ifpas/fcm and 
Pon/Poff, our analysis show that u decays exponen¬ 
tially with distance s from the free end, as shown 
in Fig. [^a)-(c). Further, the decay length of the 
exponential increases when Kp^s/km is varied from 
10-3(Figga)) to 10-2(Figgb)) and IQ-HFiggc)). 


An exponential fit to the amplitude along the com¬ 
posite for different parameter values allows us to esti¬ 
mate the decay length Ab as a function of s (Figj^d)). 
This indicates that for a fixed Pon/Poffj A is scales as 
(ATpas/^m)^- However, the decay length exhibits only 
weak dependence on motor kinetics for the range of 
attachment probabilities investigated. We also find that 
the value saturates to a constant in simulations with 
large A'pas/fcm, which we ascribe to finite system size 
effects. 

To compare with theoretically predicted decay 
length to our simulation results, we set pN = 0 in (0 
to obtain the decay length 


/ 0 ATpas \ ^ 

1 

V ^mPm J 



Since x is a function of imposed frequency and motor 
kinetics, and is to leading order independent of filament 
stiffness, the term in the square brackets is a constant 
for a given u) and Pon/Poff- As shown in FigQd), this 
prediction agrees very well with our simulations. 


Conclusions 

The scale over which active deformations persist in a 
fluctuating environment are basic questions in the study 
of living matter. Here, we address both these questions 
in the context of an ordered composite structure con¬ 
sisting of elastic elements that can be stretched and/or 
sheared by motors that can bind to and unbind from 
them, taking to account the individual motor kinetics. 
Even though extensibility is negligible locally, it affects 
the scale over which deformations persist owing to the 
competition between extensibility and shear. In partic¬ 
ular, for active systems, we show that the emergence 
of an oscillatory instability in a short segment leads to 
deformations that do not persist forever and instead 
decay over longer length scales. Further, the finite num¬ 
ber of motors attached causes an intermediate regime 
of intermittent oscillations, which was not predicted by 
models that ignore active noise. Our analysis shows that 
the feedback due to the ability of attached motors to 
sense local extension rates influences the motor activity 
and modifies the effective length scale over which strain 
decays. 

Our simulations suggest that deviations from the 
mean field theory - even for weak extensibility - can 
arise in the limit of very small number of attached 
motors (or small attachment rates) and very large num¬ 
ber of attached motors (noise acting cumulatively at 
large attachment rates). We find that, for a given motor 
activity, the decay length of strain predicted by the 
mean field theory, Ag ~ 2.{Kpas/Gpas)"^ ■, is valid even 
in the presence of noise. This confirms a finite range of 
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correlated activity that might be relevant for natural 
examples of ordered active matter such as eukaryotic 
flagella m by setting a limit over which mechanical 
coordination can persist. It also raises interesting ques¬ 
tions about how to generalize our analysis to disordered 
mixtures of motors and cytoskeletal filaments ^51 [55] i 
particularly in the limit when the coordination number 
is small and when it approaches the isostatic limit. 
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